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Single-Dirac-cone topological insulators (TI) are the first experimentally dis- 
covered class of three dimensional topologically ordered electronic systems, and 
feature robust, massless spin-helical conducting surface states that appear at any 
interface between a topological insulator and normal matter that lacks the topo- 
logical insulator ordering This topologically defined surface environment 
has been theoretically identified as a promising platform for observing a wide 
range of new physical phenomena, and possesses ideal properties for advanced 
electronics such as spin-polarized conductivity and suppressed scattering f3-22|. 
A key missing step in enabling these applications is to understand how topologi- 
cally ordered electrons respond to the interfaces and surface structures that con- 
stitute a device. Here we explore this question by using the surface deposition 
of cathode (Cu/In/Fe) and anode materials (NO2) and control of bulk doping in 
TI Bi2Se3 from P-type to N-type charge transport regimes to generate a range of 
topological insulator interface scenarios that are fundamental to device develop- 
ment. The interplay of conventional semiconductor junction physics and three 
dimensional topological electronic order is observed to generate novel junction 
behaviors that go beyond the doped-insulator paradigm of conventional semi- 
conductor devices and greatly alter the known spin-orbit interface phenomenon 
of Rashba splitting. Our measurements for the first time reveal new classes of 
diode-like configurations that can create a gap in the interface electron density 
near a topological Dirac point and systematically modify the topological surface 
state Dirac velocity, allowing far reaching control of spin-textured helical Dirac 
electrons inside the interface and creating advantages for TI superconductors as 
a Majorana fermion platform over spin-orbit semiconductors. 

The surface state of topological insulator Bi2Se3 crystals cleaved in ultra high vacuum 
(UHV) is a singly de gen erate Dirac cone that spans a bulk semiconductor band gap of 
A>300meV [sl, |6|, QJ, l20|, with a Fermi level that can be easily moved by bulk dopants 
(Fig. la). These surface electron dynamics, identified in 2008, have inspired an ongoing 
stream of proposals anticipating new phenomena and technologies that can be achieved 
by interfacing conventional electronics with two dimensional spin-chiral Dirac electrons, 
including the possibility of observing exotic quasiparticles such as Majorana fermions, dyons 
and axions with properties that have been found nowhere else in the universe [o], Q, 16 |. 
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Developments from the last few months have demonstrated that it is possi 



e to fabricate 



nanodevices using Bi2Se3 and pass current through the surface states [23|, |2J], but have 
also highlighted the fact that topological insulators are a new state of matter distinct from 
normal metals and insulators, and their properties in a junction setting are not yet known. 
Changes in the surface chemistry of Bi2Se3 have been found to modify the Dirac cone 



kinetics in unpredictable ways |2(| and induce bulk electrons to be trapped in new surface 
Dirac cones |2l| . To bridge the gap between the conceptual promise of topological insulator 
Dirac electrons and real devices that can take advantage of topological electronic order, it 
is critical to develop an experimentally motivated understanding of how interfaces between 
normal metals and topological insulators will function. 

To determine how junctions will change the topological electron kinetics, we have de- 
posited thin films on the surface of carrier doped Bi2Se3 to generate a wide range of 
semiconductor-metal P-N and N-N interfaces, and used angle resolved photoemission spec- 
troscopy (ARPES) to map changes in the topological surface electrons. The non-TI part 
of the interface is grown by evaporative deposition of strong charge donors (Cu, NO2, Fe), 
achieving conventional diode interface voltages up to Vd~0.7eV from a sub-nanometer sur- 
face layer that is permeable to photoexcited electrons [25]. The cleaved Se^~ plane of Bi2Se3 
is maintained at low temperature (T<20°K) to provide a chemically neutral and homo- 
geneous surface for these adatoms, as discussed in the online Supplementary Information 
(SI). Deposition depth is presented in units of monolayers (ML) to approximately represent 
the number of deposited ions per hexagonal unit cell of the surface, based on rough linear 
assumptions outlined in the SI. Summing the Fermi surface area of all surface states gives 
the surface state electron density or Luttinger count, which we find grows more slowly than 
the nominal rate of charge deposition with Cu"^ or Fe (Fe^"*" or Fe'^"^) ions by approximately 
one order of magnitude. After an estimated 2ML of Cu deposition the Fermi level ceases to 
move as more adatoms are deposited, and very few additional electrons can be added to the 
surface states (Fig. |2b). As shown in Fig. Ic, the topological surface state is dramatically 
changed by the voltage drop in a P-N inteface between hole doped (P-type) Bi2Se3 and 
copper, manifesting new Dirac cones and a gap in the surface density of states near the 
original (DO) Dirac point at 0.6-0.7eV binding energy. 

The observation of a gap in the surface electron density is extremely surprising, be- 
cause it is well known from theory that non-magnetic perturbations cannot open a true 
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band gap in topological surface state Dirac points . To understand how the junction 
with non-magnetic copper could cause such a feature to appear, we have created a Green's 
function implementation of the experimentally-based k-p model in Ref. [o], [ll| to numeri- 
cally simulate junctions at the surface of a semi-infinite topologically ordered Bi2Se3 slab 
based on experimentally measured bulk electron kinetics from Ref. |2ol . Unlike other 
current numerical models, this approach is capable of simulating the full screening depth 
of electrostatic interactions inside a Bi2Se3 junction (~100nm), and encompasses both the 
experimentally determined bulk kinetics and the underlying topological order of electrons 
in the simulated crystal. Using a screening potential that matches the observed change in 
topological surface state electron energies in the Bi2Se3-Cu P-N interface (A^=0.7eV), the 
simulation shows that the upper Dirac cone surface state electrons are still found in the first 
quintuple atomic layer of Bi2Se3, but the corresponding lower Dirac cone electrons are found 
more than Inm deeper in the crystal, due to hybridization with other surface states derived 
from the perturbed bulk valence band (Fig. lb). Because the lower Dirac cone electrons of 
the original surface state are deeper in the crystal than the ~lnm photoemission penetration 
depth, photoemission measurements are expected to show a void or "gap" near the DO Dirac 
point, even though the three dimensional band structure is not gapped. This effect observed 
from our data and model is consistent with a recent theoretical proposal that interfacing TI 
surface states with a parasitic ( deg enerate) 3D metal can reshape the surface bands to take 
the appearance of a band gap [27], a coexistence expected in many theoretically predicted 
topologically ordered metals such as half-Heusler compounds [28 1. 

The incremental evolution leading to the appearance of multiple new Dirac cones on 
Cu-deposited P-type Bi2Se3 is shown in Fig. [2l The build-up of positively charged surface 
ions first donates electrons to the original Dirac cone of freshly cleaved Bi2Se3, and then 
gradually binds four new surface states that intersect in two new Dirac points (Dl and D2) 
at the Brillouin zone (BZ) center. Tracing these new states from their points of intersection 
in the BZ center, we observe that they exchange partners and merge with a different surface 
band at momenta larger than ~0.lA~^, corresponding roughly to the outer contour of the 
region in which symmetries are inverted to bring about the TI state in numerical models. 
The resulting Fermi surface after heavy deposition (Fig. [2^, right) is made up of three 
rings, of which the outer two are doubly degenerate (within intrinsic width resolution). 
The doubly degenerate outer rings cannot carry a non-zero Berry's phase, meaning that 
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the Berry's phase of n by which a topological insulator interface is identified is obtained 
from the innermost singly degenerate Fermi surface ring, located deepest inside the Bi2Se3 
crystal. The outermost ring is found at k~1.7-2.0A~^, and is hexagonally warped into a 
doubly-degenerate "snowfiake" shape that resembles the singly degenerate hexagonal Fermi 
surface associated with magnetic instabilities in related Bi2Te3 [llj]. Although alternative 
terminologies have been proposed in recent literature 29|], we refer to all of these surface 
bands as topological surface (or interface) states, because they manifest at the surface of 
the material and all five Dirac bands must be observed to measure the 0B=7r topological 
Berry's phase. 

In topologically trivial semiconductors with strong spin-orbit interactions, electrons near 
an interface are commonly spin-split into spin polarized Dirac cones by the Rashba effect. 
The partner swapping behavior we observe in the Dl and D2 Bi2Se3-Cu P-N interface 
states appears identical to the conventional Rashba effect in photoemission measurements 
at small momenta (K^O.lA"^), but has very different phenomenological imphcations. The 
conventional Rashba effect causes surface electrons to be more strongly spin-polarized as 
their momentum deviates from the BZ center, due to the Rashba Hamiltonian term which 
is proportional to momentum {Hji={2aRk x z) ■ s, with s representing spin). In contrast, 
though Bi2Se3-Cu P-N interface electrons obey this linear relation near their Dirac points, at 
large momentum the surface electrons are all nearly doubly-degenerate, and do not generate 
a large density of spin-polarized conduction states at the Fermi level. The difference can be 
seen clearly by comparing the Bi2Se3-Cu interface band structure with the surface states of 
BiTel (Fig. |2]:), a compound chemically similar to Bi2Se3, but which our measurements show 
to be topologically trivial and have an extremely strong Rashba effect (ai?=2.2±0.1eV*A, 
see SI). Due to the partner swapping effect, spin-helical conduction states inside the interface 
are still limited to a single Dirac cone as is desired for applications. 

The upper Dirac cone surface state of the N-type Bi2Se3 samples explored in this study 
(including Cu^Bi2Se3) has a momentum-axis width {6k) of <5k~0.035A~^ at half maximum 
intensity before deposition, which can be inverted to reveal the electronic mean free path 
(mfp) of i=l / 6k=29A. Momentum- axis width of the original Dirac cone increases steadily 
as Cu is added to the surface (Fig. [2li) and is surprisingly broad when one considers 
that disorder-induced scattering is suppressed by the spin-momentum locking of electrons 
0, y, [lol. Depositing magnetic iron ions instead of non- magnetic copper results in an 
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unusual fluctuation in 5k that has been discussed in Ref. [2l| and may be due to reduced 
disorder following a 2D ferromagnetic ordering transition. The new surface states generated 
by copper deposition have sharper momentum widths of ~0.027±3A^^ {£=37A mfp) just 
above the 0.27eV binding energy Dl Dirac point in Fig. |2^(right) and ~0.023±3A"^ {i=43A 
mfp) just above the O.leV D2 Dirac point. The capability to generate a longer mean free 
path of spin polarized surface electrons via deposition or top gating, as identified in Fig. |2l 
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22| and non-Abelian coherent behaviors 



is desirable to enhance long range rnagnetic 
unique to topological insulators jo], Q . 

Topological insulators are among the very few material systems with band structures that 



can be manipulated to achieve transport defined by a Dirac point [5[, and most proposed 
applications for TIs derive from the behavior of low density spin-helical Dirac surface elec- 
trons. In Fig. [St we present two realizations of P-N interfaces in which the Bi2Se3 surface 
state Dirac point is brought to the Fermi level where it will shape transport properties. The 
surface chemical potential of the N-type sample shown at left (as-grown Bi2Se3 with [Se]^"*" 
defects) is lowered by depositing N02~, and the P-type sample (Bii.ggsCao.oosSea) was mod- 
ified through deposition of In'*". Our data show that lowering the surface chemical potential 
of an N-type sample to the Dirac point has the effect of dramatically increasing the Dirac 
velocity from V£i=1.55±0.08eV-A to 2.30±0.08eV-A, observed from the lower Dirac cone 
band slope traced in Fig. [3^,c. The velocity of lower Dirac cone electrons changes rapidly 
immediately after deposition begins, but stabilizes after the surface chemical potential falls 
below the bulk conduction band minimum 0.17±0.02 eV above the Dirac point (Fig. [3b). 
When the chemical potential lies inside the bulk electronic gap as measured by surface- 
sensitive photoemission, the lower Dirac cone velocity has the same value of 2.3±0.1eV-A 
for both N-type and P-type samples (Fig. [3)d,c) even though the they are exposed to surface 
ions with opposite charge. As the chemical potential at the surface of Bi2Se3 is bent into 
the bulk band gap, numerical simulations show that the region near the surface becomes 
depleted of bulk conduction electrons (Fig. [HJi and Fig. S5 in the SI), causing bulk and 
surface charge carriers to no longer spatially overlap with one another. The bulk and surface 
charge carriers occupy the same atomic orbitals, and spatial overlap between the surface and 
bulk conduction electrons will result in strong on-site charge and spin interactions that are 
likely responsible for the tunable Dirac velocity observed in our data. 

Interfaces with like carrier type (N-N or P-P) are also of great technological importance. 



6 



because they create a less resistive electronic junction with minimal charge inhomogeneity. 
Measurements in Fig. H] explore the effect of light copper deposition on the surface of N-type 
Cuo.o2Bi2Se3 to observe how an N-type surface layer interfaces with the electron doped topo- 
logical insulator. Excess surface-deposited Cu adatoms cause the electron energies near the 
surface of the TI to bend towards higher binding energy, as they would in an N-N Schottky 
diode. A new surface state is fully formed after 0.08ML of Cu deposition (Fig. H^, mea- 
surements at <0.08ML coverage in SI), and remains doubly degenerate within experimental 
resolution until approximately 0.15ML of Cu deposition, at which point it becomes spin-split 
as seen from a deviation between the band energy at the Brillouin zone center and the band 
energy minimum (Fig. Ht). 

Because TI surface states cross the bulk band gap and have different energetics than 
bulk electrons ^ 0], when interfacing non-TI materials with a topological insulator, the 
surface states will trap charge carriers and attract or repel bulk electrons even when the 
bulk carrier doping of each material is identical. Our ARPES measurements show that 
prior to depositing ions on the vacuum-cleaved Cuo.o2Bi2Se3 surface, a negative charge is 
trapped by the surface state, and will be screened over a range of more than 10 nanometers 
as simulated from the Poisson equation (Fig. U^). Assuming that each Cu atom enters a 
+1 valence state, depositing 0.08ML of copper causes the net charge at the crystal surface 
to change from negative to positive. Between and 0.08ML (shaded in Fig. H^) the electric 
field penetrating beyond the first few nanometers of the crystal surface is expected to be 
very small, allowing electrons to retain their bulk-derived energetics and symmetries as 
they approach the surface from deeper in the crystal. As Cu deposition is raised to 0.16ML, 
electron density in the surface state grows more slowly than the rate of deposition, causing an 
increasing proportion of the surface charge to be screened by bulk electrons. Titrating with 
positive ions or a gating voltage (~0.1eV from our data) to achieve an electrically neutral 
surface is therefore a promising method to enhance bulk interactions with the surface states, 
such as magnetic and superconducting proximity effects. 

Recent measurements have identified that the bulk doped TI Cua;Bi2Se3 (x>0.1) is a 
type-2 N-type bulk superconductor (Ref. |20| and Fig. la, left), and that magnetic vortices 
at the surface are likely to bind non-Abelian particles known as Majorana fermions that have 
not yet been observed in nature jo], [2^ . Our calculations in Fig. show that an electric field 
from the negatively charged surface state will cause rapid phase fluctuations that suppress 
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superconductivity if the surface charge is not neutrahzed (Fig. H^, inset). Moreover, because 
vortices are repelled by regions of high superfluid density, this surface environment provides 
a new mechanism to control the dynamics of Majorana fermions and potentially observe 
their non-Abelian properties. For normal superconductors, surface perturbations suppress 



superfluid density, generating an attractive force that pins vortices in place [30|, |3l| • For a 
topological insulator however, our data show that adjusting the surface chemical energy to 
cancel out the negative surface state charge can enhance superfluid density under a probe, 
causing surface vortices to be repelled by the probe rather than pinned to it (Fig. HJi). 
The control and exclusion of vortices is a subject of active research to improve the proper- 



ties of superconductors [17|, [SiJ, and the intrinsically charged TI surface provides a new and 
promising platform for vortex manipulation. Taken collectively, the P-N and N-N interface 
electronic structures observed here reveal how interface effects such as Rashba spin-splitting 
and the superconducting proximity effect are changed by the presence of topological order, 
and how topological insulator interfaces can be tuned to adjust the average speed of surface 
Dirac electrons and interact with surface Majorana fermion vortices. 

Methods summary: Angle resolved photoemission spectroscopy (ARPES) measure- 
ments were performed at the Advanced Light Source beamlines 10 and 12 using 35.5-48 
eV photons with better than 15 meV energy resolution and overall angular resolution bet- 
ter than 1% of the Brillouin zone (BZ). Samples were cleaved and measured at 15°K, in a 
vacuum maintained below 8x10^^^ Torr. Fe atoms were deposited using an e-beam heated 
evaporator at a rate of approximately O.lA/minute, and Cu atoms were deposited at a rate 
of 0.4A/minute. A quartz micro-balance supplied by Leybold-Inficon with sub- Angstrom 
sensitivity was used to calibrate the deposition flow rate, and a factor of 2.5A/ML was used 
to estimate monolayer depth. Surface charge is estimated from the expected -|-1 ioniza- 
tion state of copper on the cleaved [111] Se surface [2^. Adsorption of NO2 molecules on 
Bi2Se3 was achieved by controlled in situ exposures under static flow mode, with care to 
minimize photon exposure of the adsorbed surface. Crystal growth of large single crystals 
of Ca2:Bi2-xSe3 and CuxBi2Se3 and numerical techniques used to simulate the perturbed TI 
surface are described in the SI. 
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FIG. 1: Coulomb-biased topological surface states: a, The (blue) bulk and (green) surface 
band structure of Bi2Se3 is traced above ARPES measurements of band structure at n-type and 
p-type bulk compositions, b, A numerical simulation of copper-interfaced Bi2Se3 shows the partial 
density of low energy surface states at depths of 0-0. 5nm and 2-3nm inside the crystal, (white 
dashed lines) Electrons from the upper Dirac cone of 612863 remain in the outermost 0.5nm, and 
electrons in the lower Dirac cone are located l-2nm deeper in the crystal, c, The surface states 
of Cu-deposited Bi2Se3 resonate at different photon energies, demonstrating that they occupy a 
quantum well-like environment. Measurements have a ^Inm penetration depth into the crystal 
25l |. An integrated comparison reveals the band structure traced at right, with new Dirac points 
labeled Dl and D2 and the Fermi surface shown as an inset. States that evolve adiabatically from 
the original upper Dirac cone ("DO") are drawn in a darker hue. 
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FIG. 2: A topological p-n interface: a, Surface deposition of copper is shown on bulk p- 
type Bi2Se3 (0.5% Ca doped as Bii.ggsCao.oosSes), with approximate pre-deposition electronic 
structure traced on the first panel and energies shifted to align surface band features between 
panels. The mean free path estimated from inverse momentum width is labeled in nanometers for 
surface conduction bands. Fermi surfaces are shown for (inset) 0.25ML and (right) 2.2ML of Cu 
deposition, with numerically predicted spin orientations labeled with arrows, b, The surface state 
carrier density estimated from photoemission data (Luttinger count) is plotted as a function of Cu 
deposition, c, (top) Large Rashba splitting in the surface states of topologically trivial BiTel is 
compared with (bottom) the strikingly different Rashba-like effect in topologically ordered surface 
electrons, d, The imaginary self energy of the original ('DO') upper Dirac cone after Cu and Fe 
deposition is obtained by multiplying the momentum-axis width {5k) by half of the surface state 
velocity (Im(S)=^^^^^^f^). Deposition depth axes are mismatched in proportion to the expected 
surface valence states of Cu"'" and Fe'^+. 
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FIG. 3: Dirac point P-N interfaces: a, A measurement of band structure in N-type as-grown 
Bi2Se3 is labeled with the lower Dirac cone velocity and Dirac point energy (fJ-o) before surface 
deposition of NO2. (b) The lower Dirac cone band velocity is plotted as a function of chemical 
potential above the Dirac point as NO2 is added to the surface of N-type Bi2Se3. c, Photoemission 
images of bulk P-type and N-type samples with chemical potential set at the Dirac point via 
surface deposition of negatively and positively charged molecules, respectively (NO2 and In). An 
inset shows the Dirac node Fermi surface, and the lower Dirac cone velocities (slopes) of the pristine 
pre-deposition surface states are traced in white, d. The electron density distribution is modeled 
for (gray) as-grown Bi2Se3 and (gold) the same sample with chemical potential set at the Dirac 
point by NO2 deposition. 
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FIG. 4: Superconductors and n-n junctions: a, Photoemission measurements show the evo- 
lution of surface electron kinetics with light Cu surface deposition on bulk-doped Cua;Bi2Se3. b, 
The new band that appears after 0.08 ML Cu deposition, labeled 'Dl' at the BZ center, has no 
momentum dispersion along the z-axis. The original surface state Dirac point is labeled 'DO', c. 
The doping evolution of band energies at the T-point is shown. The Dl band minimum, labeled 
'Dl-min', and Dirac point energy ('Dl') diverge after heavy deposition, d. Probes interfacing 
with normal superconductors suppress the superfluid density, attracting ("pinning") vortices. On 
a topological insulator, a gentle probe with bias voltage similar to O.leV can attract vortices by 
enhancing the superfluid density. Regions with suppressed superfluid density are indicated with 
green shading, e. The z-axis electron density distribution is modeled for different degrees of surface 
doping, (shaded region) Deposition levels between OML and 0.08 ML will reduce bulk screening, 
strengthening superconductivity near the surface. An inset shows the electric field inside the crys- 
tal. 
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